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Abstract. In the gravitational lens system B1600-)-434 the brighter image. A, is known to show rapid variability 
which is not detected in the weaker image, B (Koopmans & de Bruyn 2000). Since correlated variability is one 
of the fundamental properties of gravitational lensing, it has been proposed that image A is microlensed by stars 
in the halo of the lensing galaxy (Koopmans & de Bruyn 2000). We present VLBA observations of B1600-I-434 
at 15 GHz with a resolution of 0.5 milliarcsec to determine the source structure at high spatial resolution. The 
surface brightness of the images are significantly different, with image A being more compact. This is in apparent 
contradiction with the required property of gravitational lensing that surface brightness be preserved. Our results 
suggest that both the lensed images may show two-sided elongation at this resolution, a morphology which does 
not necessarily favour superluminal motion. Instead these data may suggest that image B is scatter-broadened at 
the lens so that its size is larger than that of A, and hence scintillates less than image A. 
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1. Introduction 

The gravitational lens B1600-f434 was discovered by 
Jackson et al. (1995). The radio source is a quasar at 
z = 1.59 and has two images separated by 1.4 arcsec, 
which are denoted A and B. The main lensing galaxy 
is identified as an edge-on spiral at z = 0.41, with an 
additional lensing contribution from a nearby companion 
galaxy (Jaunsen & Hjorth 1997, Fassnacht & Cohen 1998, 
Koopmans, de Bruyn & Jackson 1998). The background 
quasar is variable in both optical and radio bands, thus 
allowing time delay measurements from monitoring obser- 
vations (Jaunsen & Hjorth 1997, Koopmans et al. 2000). 
A time delay of 47lg^ days has been determined from the 
radio monitoring (Koopmans et al. 2000), and a delay of 
51 ± 4 days from the optical observations (Burud et al. 
2000). 

The radio observations of Koopmans et al. (2000) 
showed that the brighter image. A, varied more rapidly 
than the weaker image, B. The rapid variability in A was 
not found to be correlated with corresponding rapid vari- 
ability in B, although differential variability between A 
and B on longer time scales allowed the measurement 
of the time delay mentioned above. Correlated variability 
is one of the fundamental characteristics of gravitational 
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lensing. Additional variability in A can only be explained 
by invoking propagation effects which affect image A only 
(Koopmans & de Bruyn 2000; hereinafter KdB). Such ef- 
fects include a combination of scattering in the lensing 
galaxy and scintillation in our own galaxy, or microlens- 
ing in the halo of the lensing galaxy. KdB considered the 
possibility that B might be scatter-broadened in the lens 
and thereby would scintillate less than A. They argued 
however that the time scale and the frequency dependence 
of the variability, taken in conjunction with the high scat- 
tering measure required in the lens galaxy, did not favour 
scintillation as a possible cause. Instead, they proposed 
microlensing of image A in the halo of the lensing galaxy 
as a preferred alternative if the radio source has a super- 
luminal component. 

One of the ways to test these scenarios is to image 
the structure in the lensed images at high spatial reso- 
lution, and to apply the constraints imposed by gravita- 
tional lensing. In particular, the surface brightness of grav- 
itationally lensed images must be equal and any significant 
departure would suggest that one of the images has been 
affected by scattering more than the other. Since the mi- 
crolensing hypothesis requires the background quasar be 
a superluminal source, the images, therefore, might be ex- 
pected to show evidence of this motion in their morphol- 
ogy, such as a core-jet structure typical of such quasars. 
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In addition, the polarization properties of the images can 
differ if one of the ray paths is affected either by scattering 
or microlensing. 

With the above aims in mind, we observed B1600+434 
at 15 GHz using the VLBA with a resolution of 0.5 mil- 
liarcsec (mas), to deteminc the fine-scale structure and to 
measure image sizes for A and B. The observations are 
described in Section 2 and the results are discussed in 
Section 3. 



2. Observations and data analysis 

We observed B1600+434 on 2000 September 9/10 (1900 
UT to 0500 UT) at 15 GHz using all ten telescopes of 
the VLBA. The data were recorded in eight 8-MHz chan- 
nels in dual circular polarization, yielding a total recorded 
bandwidth of 32-MHz per polarization. Two-bit sampling 
quantization was used. A single correlation centre was se- 
lected at a = 16^01"40.47^ 6 = 43'^16™47.2^ (J2000). 
Since the image separation of 1.4 arcscc is large compared 
with the resolution of 0.5 mas, it is necessary to sample the 
data finely both in frequency and time to avoid smearing 
of the visibility function and the associated image distor- 
tions in the field of view. The data were therefore sampled 
at 0.5 MHz in frequency and 0.5 sec in time. 

We observed 3C345 (J1640+3946), 1611-1-343 
(J1613-F3412) and 1749-^096 (J1751-F0939) as cali- 
brators. The total integration time on the target source 
B1600+434 was 195 min. Wc used 1749+096 to calibrate 
the instrumental polarization, as well as to determine the 
bandpass function. 

The method of data analysis follows closely that de- 
scribed by Patnaik et al. (1999) as used for the gravi- 
tational lens B1422-I-231. The data were analysed using 
the software package AIPS, maintained by the National 
Radio Astronomy Observatory (NRAO). After correcting 
for the parallactic angle variation of the telescopes, am- 
plitudes were calibrated using the system temperatures 
and telescope gain factors as a function of elevation. The 
phase slopes and offsets across each 8-MHz band were 
aligned using the pulsc-cal information from each tele- 
scope. Standard fringe-fitting was performed on the cali- 
brators. After careful editing, the bandpass function was 
determined using the source 1749+096. Wc followed the 
usual procedure for calibrating instrumental polarization, 
using 1749+096 as the polarization calibrator. Typical in- 
strumental polarization amplitude was measured to be 
~1.6%. Since B1600+434 has two images with similar fiux 
densities at a large angular separation, we had to iterate 
the fringe-fitting process several times, each time using 
the model of the source determined from the previous so- 
lution. 

The final map was made with two sub-fields, centred 
on each lensed image. We also made polarization maps of 
the source in Stokes /, Q and U. 



Image Peak Total 
in mjy/beam in mjy 



A 
B 



9.0+0.15 
6.6+0.15 



17.2+0.5 
13.5+0.5 



Size 

mas X mas, PA 

0.69 X 0.25, 103° 
0.73 X 0.30, 93° 



Table 1. Peak (in mJy/beam) and total (in mJy) fiux 
densities of the lensed images of B1600+434 at 15 GHz. 
The image deconvolved sizes were obtained by fitting a 
single gaussian in the program JMFIT. The error in the 
measured size is 0.01 mas in each axis for image A and 
about 0.015 mas for image B, and 1° in position angle. 
The image separation of B from A is 723.08+0.05 mas to 
the east (in RA) and 1189.16+0.05 mas to the south (in 
Dec). 



3. Results and Discussion 

The maps of the lensed images, made using uniform 
weighting, are shown in Fig.l. The maps have been re- 
stored using a circular gaussian with a FWHM of 0.5 mas. 
We did not detect any polarized emission from either com- 
ponent A or B, yielding a 5a limit to the percentage linear 
polarization of 2.6% for image A, and 3.3% for image B. 
The rms noise in the polarization map is 0.09 mJy/beam. 

Table 1 lists the results of fitting a single gaussian com- 
ponent to each of the images. Two clear results emerge 
from these observations. First, both the images are re- 
solved, and A is more compact than B. Secondly, the struc- 
ture is elongated on cither side of the peak in both of the 
images. In order to quantify the errors of the size mea- 
surements we divided the data into two sets by removing 
alternate scans so as to preserve similar uv-covcragc. The 
resulting maps were then fitted with a single gaussian. The 
results show that the errors for image A are ~ 0.01 mas 
in each axis. However, the errors arc larger for image B 
by about 50%. Note that the error estimates derived in 
this way fold in most systematic calibration or imaging 
errors which may affect the component size measurement, 
excluding deviation of the component from non-Gaussian 
shape. A weak extension is seen to the south of image B 
which appears to be genuine as it is present in both the 
datasets mentioned earlier. We are not entirely certain if 
it is real as there is no counterpart in image A. This can be 
confirmed by another independent dataset. The measured 
size of image B, however, is not affected by including this 
feature. 

The implication of the first result is that the surface 
brightnesses of the two images differ significantly. The fiux 
density ratio of A/B is 1.27+0.05 whereas the ratio of im- 
age sizes is 0.80+0.15. Secondly, the image morphology as 
measured is in apparent conflict with a typical superlumi- 
nal radio source in the sense that such radio sources are 
core-dominated and must have one-sided jets due to rela- 
tivistic motion Doppler boosting of the knot components. 
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Fig. 1. 15 GHz map of B1600+434 A (left) and B (right). The contour levels for both the images are -2, -1, 1, 2, 
4, 8, 16, 32, 64, 128 x the contour interval which is chosen as the Scr noise in the map. For image A, contour interval 
is 0.42 mJy beam~^ and peak flux density is 9.5 mJy beam~^ and for image B contour interval is 0.42 mJy beam~^ 
and the peak is 6.7 mJy beam~^. The convolving beam of 0.5 mas circular gaussian is drawn at the lower left-hand 
corner in each map. 



In order to understand the higher modulation index 
(fractional rms variability) of image A (2.8% at 8.4 GHz), 
KdB considered the possibility that the size of B (modu- 
lation index of 1.6% at 8.4 GHz) might be larger through 
scatter-broadening by the ionised interstellar medium in 
the lens galaxy. The fact that the surface brightnesses 
do not agree, and in particular, that the weaker image 
B is more extended, suggests that B has been affected 
by propagation through the lensing galaxy more than A. 
This would support the idea that B is somewhat broad- 
ened by scattering in the lensing galaxy. We note that its 
ray path passes close to the bulge of this edge-on spiral 
galaxy whereas the ray path of A goes through the halo of 
the lens. For the two images to scintillate in our Galaxy, 
the required sizes are 62 /uarcsec and 108 /uarcsec for A 
and B respectively (KdB). Our observations cannot re- 
solve such sizes and our measured sizes are considerably 
larger than these. We note, however, that the observed 
images have been parametrised as single gaussian com- 
ponents. This does not exclude any sub-structure in the 
source. Therefore, our result that image B is larger in size 
than image A is in general agreement with the expectation 
from scintillation taking place in our Galaxy. 

KdB measured a modulation index <1.2% at 1.4 GHz 
and 3.7% at 5 GHz. If the variability was due to scin- 
tillation, then the modulation index should have in- 
creased going to lower frequencies. However, if there is 
frequency-dependent structure then the modulation index 
can be smaller at lower frequencies. A flat spectrum source 
likely consists of more than one self-absorbed component 
(Cotton et al. 1980). It is then possible that the compo- 



nent which is compact and scintillates at 5.0 GHz, which 
may for example be the core, does not contribute much to 
the flux density at 1.4 GHz due to self-absorption. Some 
other component in the jet might be brighter at 1.4 GHz 
but this component would be more extended and therefore 
unlikely to scintillate. Therefore, a smaller modulation in- 
dex at a lower frequency does not necessarily rule out scin- 
tillation occurring at a higher frequency. Such frequency- 
dependent structure has been noted in the gravitational 
lens B0218-h357 (Patnaik & Porcas 1999). 

The alternative explanation proposed by KdB is that 
A is microlensed in the halo of the lensing galaxy. However, 
for the microlensing to match the time-scale of variabil- 
ity, the model used by KdB requires that the background 
source be a supcrluminal source with a knot the size of a 
few yuarcsec containing ~ 5-11% of the flux density, moving 
at an apparent speed in the range ~ 9-26c. The morphol- 
ogy of a supcrluminal source will of necessity be that of 
a bright core and one-sided jet, due to Doppler boosting 
of a relativistically moving knots in the jet. The detec- 
tion of extensions on either side of the peak in both the 
lensed images here is contrary to the above expectation. 
It would, therefore, appear that supcrluminal motion in 
the background source is not favoured and hence the mi- 
crolensing caustic crossing time will be much longer than 
the observed variability. However, it is also possible that 
the peak we observe is not necessarily the core; it might be 
a bright knot and thus still allow the possibility that the 
source could be intrinsically one-sided. We do not have 
sufficient resolution to distinguish this level of detailed 
structure. If we assume that the peak corresponds to the 
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core, which may be a reasonable assumption given the fact 
that these observations are at 15 GHz and with 0.5 mas 
resolution where cores usually brighter than knots, then it 
would impliy that microlensing combined with superlumi- 
nal motion in the background source may not be a viable 
explanation for the rapid variability of the lensed image 
A. 

As an example, the two brightest images in the lens 
system B1422+231 have elongation produced by lensing 
(Patnaik & Porcas 1998, Patnaik et al. 1999). In this lens 
system the background radio source is also a flat spectrum 
object and at milliarcsec-scales does not have the standard 
core-jet structure. It has been suggested that the magni- 
fication in this source might be very high and that the 
background object could possibly be a radio-quiet object. 
A similar scenario is possible for B1600-I-434. Certainly, 
in both B1422-I-231 and B1600+434, no jet has been de- 
tected and both sources show two-sided emission. If the 
extension of the images we detect is due to stretching of 
the extended background emission by gravitational lens- 
ing, then we cannot be sure of the direction of the jet. 
However, should there be superluminal motion (which can 
be as much as 0.15 milliarcsec per year for the bright 
quasars), we would detect changes in image size after a 
year or two. 
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4. Conclusions 

We have observed the gravitational lens system 
B1600+434 at 15 GHz using the VLBA at a resolution of 
0.5 milliarcsec. We find that the two lensed images are re- 
solved and have elongation on either side of the peak, and 
that the surface brightnesses of the images arc significantly 
different. The image shapes do not appear to favour super- 
luminal motion in the background radio source and hence 
the interpretation that A is microlcnscd is not favoured by 
these observations at this epoch. Our results are consistent 
with the view that image B is perhaps scatter-broadened 
by the interstellar medium of the lens galaxy, and there- 
fore scintillates less than image A in our Galaxy. However, 
VLBI observations of the source at a range of other ob- 
serving frequencies are required to provide further insight 
into this question. 
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